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We propose a scheme for laser cooling of negatively charged molecules. We briefly summarise the 
requirements for such laser cooling and we identify a number of potential candidates. A detailed 
computation study with Cf, the most studied molecular anion, is carried out. Simulations of 3D 
laser cooling in a gas phase show that this molecule could be cooled down to below 1 mK in only 
a few tens of milliseconds, using standard lasers. Sisyphus cooling, where no photo-detachment 
process is present, as well as Doppler laser cooling of trapped C^, are also simulated. This cooling 
scheme has an impact on the study of cold molecules, molecular anions, charged particle sources 
and antimatter physics. 

PACS numbers: 37.10.Mn, 37.10.Rs 


Molecular anions play a central role in a wide range 
of fields: from the chemistry of highly correlated sys¬ 
tems w to atmospheric science to the study of the 
interstellar-medium nu. However, it is currently very 
difficult to investigate negative ions in a controlled man¬ 
ner at the ultracold temperatures relevant for the pro¬ 
cesses in which they are involved. Indeed, at best, tem¬ 
peratures of a few kelvins have been achieved using su¬ 
personic beam expansion methods or trapped particles 
followed by electron cooling, buffer gas cooling or resis¬ 
tive cooling hqhhi. The ability to cool molecular an¬ 
ions to sub-K temperatures would finally allow investiga¬ 
tion of their chemical and physical properties at energies 
appropriate to their interactions. Furthermore, anionic 
molecules at mK temperatures can also play an impor¬ 
tant role in (anti)atomic physics, where copious produc¬ 
tion of sub-K antihydrogen atoms currently represents 
the dominant challenge in the field. Sympathetic cool¬ 
ing of antiprotons via laser-cooled atomic negative ions 
that are simultaneously confined in the same ion trap has 
been proposed as a method to obtain sub-K antiprotons 
15] . As an alternative to this yet-to-be-realized proce¬ 
dure, ultra-cold molecular anions could replace atomic 
anions in this scheme, and would thus facilitate the for¬ 
mation of ultra-cold antihydrogen atoms. More generally, 
cooling even a single anion species would be the missing 
tool to cool any other negatively charged particles (elec¬ 
trons, antiprotons, anions) via sympathetic cooling. 

In this letter we present a realistic scheme for laser 
cooling of molecular anions to mK temperatures. Laser 
cooling of molecules has been achieved only for a very 
few neutral diatomic molecules (SrF, YO, CaF) p~6l - H8] . 
Furthermore, even if well established for neutral atoms 
and atomic cations, laser cooling techniques have so far 
never been applied to anions [T9j. This is because in 
atomic anions, the excess electron is only weakly bound 
by quantum-mechanical correlation effects. As a result, 
only a few atomic anions are known to exhibit electric- 
dipole transitions between bound states: Os - , La - and 
Ce - [IM2]. 

For molecules the situation is quite different because 


their electric dipole can bind an extra electron, and even 
di-anions have been found to exist [23]. For instance, 
polar molecules with a dipole exceeding 2.5 Debye exhibit 
dipole-bound states. Highly dipolar molecules such as 
LiH - , NaF - or MgO - possess several such dipole-bound 
states [24] US], and valence anionic states exist as well. 
For simplicity, in this letter we only focus on diatomic 
molecules, even if the Sisyphus laser cooling techniques 
proposed here can also be applied to trapped polyatomic 
molecules [26] . In the Supplementary material, Table [I] 
[27] . we present a review of most of the experimental as 
well as theoretical studies of diatomic anions, with useful 
references, if further studies are required. 

The first excited state (and sometimes even the ground 
state) of many molecular anions he above their neutral¬ 
ization threshold, such as in the case of H</, CO - , NO - , 
Ng , CN - or most of systems with 3, 4 or 11 outermost 
electrons. For this reason, these anions are not stable 
against auto-detachment processes and exhibit pure ro- 
vibrational transitions with ^ 100 ms lifetimes. Even if 
such long-lived states can still be of interest for Sisy¬ 
phus cooling, for narrow-line cooling or for Doppler laser 
cooling in traps £6 [28] , pure electronic transitions are 
preferred for rapid laser cooling. Transitions of ^ 100 ns 
lifetime can be found between well-separated electronic 
states (typically Be^X states), whereas transitions in 
the infrared region between electronic states (typically 
Ae^X states) have longer lifetimes of ~ 100 ps. Some 
challenges in laser-cooling molecular ions have been de¬ 
scribed in [29] . For direct laser cooling of neutral di¬ 
atomic molecules, a key ingredient is a good branching 
ratio, i.e. Franck-Condon factor (FC factor), between 
vibrational levels (SrF, YO, CaF all have more than a 
98% branching ratio on the A 2 n 1 / 2 (^ / = 0, J r = 1/2) 
X 2 L(v" = 0,7V 7 ' = 1) transition). Even with these con¬ 
siderations, the choice of the most suitable candidate is 
not obvious and is a compromise between using fast elec¬ 
tronic transitions and choosing extremely good FC fac¬ 
tors. Indeed, molecules with good FC factors can be 
found for weak dipole-valence bound transitions or for an¬ 
ionic molecules with 6 or 12 outermost electrons but with 
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forbidden dipole transitions. FC factors greater than 70% 
can be found in systems with 8, 14 or higher numbers of 
outermost electrons or in molecules which include Li or A1 
atoms [SO] [31]. Unfortunately the corresponding transi¬ 
tions are often in the deep infrared region. The best com¬ 
promise seems to be the systems with 9 outermost elec¬ 
trons having X(cr 1 7r 4 ) 2 l^, A(cr 2 7t 3 ) 2 n( n ), Bfy 2 ?! 4 ) 2 !^ 
levels. A list of such systems is given in Table [l| 

Clearly the most studied system is C^~, with a perfectly 
known spectrum. exhibits a B 2 L(v' = 0) <n>X 2 L(v" = 
0) system and a A 2 n 1 / 2 (^ / = 0) X 2 L(v" = 0) system 
with 72% and 96 % branching ratio, respectively. Be¬ 
sides, this anion has the notable advantage of not pre¬ 
senting any hyperfine structure. As a potential further 
benefit of this system, we mention that through laser 
photo-detachment of cold C^", we could produce cold 
C 2 molecules, important in combustion physics and as¬ 
trophysics. Even if more studies on laser cooling are 
needed, C 2 looks like a suitable candidate to be further 
laser cooled near 240 nm on the 0-0 Mulliken band (d 
1 L+ X 1 ! 4- ) which has an extremely favorable branch¬ 

ing ratio of 99.7 % [32] . 


| group | example 


I-VII 

LiF- LiCl- NaF“ NaCl" MnH-(X 6 A) 

II-VI 

ZnO - BeO - MgO - ZnF - 

III-V 

BP- AIN- A1P- AlAs- GaP- InP“ GaAs - BN“ GaN - 

IV-IV 

C 2 Si^ CSi-(B unstable) Sn“ Pb^ SnPb“(X 2 n) 


TABLE I: Diatomic anions with 9 outermost electrons, listed 
in terms of group of atoms. A more complete table with all 
diatomic molecules, and references, is given in the Supplemen¬ 
tary material. They are isoelectronic to the neutral molecules 
CN and SrF with the same X 2 1 + A 2 TT B 2 L + electronic struc¬ 
ture with sometimes stable quartet excited states. 

We will therefore concentrate on C</ as a benchmark 
to study laser cooling of anionic molecules. Note however 
that several other molecules, such as BN - or AIN - , may 
be used as well. They offer very similar structure prob¬ 
abilities with better FC factors (higher than 98%) but 
with a B 2 I —>A 2 n decay channel that is absent in the 
homonuclear case of C^~. Contrary to C^~, heteronuclear 
molecules present a closed rotational level scheme. How¬ 
ever, further spectroscopic studies are clearly required for 
such systems, as well as for other promising ones, such 
as metal-oxyde systems (FeO - , NiO - ) or hydride ones 
(CoH - ). 

In order to study laser cooling of C^~ we have performed 
three types of simulations. The first one is ” standard” 
laser cooling in the gas phase (thus with no strong ex¬ 
ternal fields present); the second one is laser cooling in a 
Paul trap; and the last one is Sisyphus cooling of trapped 
ions in a Penning trap. The simulations are performed 
with the C++ code described in [28], which now also in¬ 
cludes full N-body space charge effects [33]; the Lorentz 
force is solved by using Boris-Buneman integration algo¬ 
rithms om>]. Briefly, a Kinetic Monte Carlo algorithm 


gives the exact time of events (absorption or emission 
of light) when solving the rate equations to study laser 
excitation of the molecules under the effect of Coulomb, 
light scattering, dipolar, Stark and Zeeman forces. The 
C^ energy levels and required laser transitions are shown 
in Fig. [^a). 



FIG. 1: (Colour online). C^ level structure and laser cool¬ 
ing transitions, (a): Electronic and vibrational levels includ¬ 
ing the photo-detachment threshold of neutral C 2 ( 1 Lg , 3 IT U ), 
from m- The widths of the arrows are proportional to the 
transition strengths given by the FC. The green dotted line 
indicates that photo-detachement can occur from the B state 
due to a second photon absorption, (b): Detail of the Sisy¬ 
phus cooling principle with a zoom on the X o A energies: 
|Xi}=X(i/' = 0, N" = 0) and |X 2 )=X(u" = 0,1V" = 2), 
|A)=A(+ = 0, N' = 1, J' — 1/2). A Penning-like trap is pre¬ 
sented in the upper part. A magnetic gradient field creates 
the Sisyphus potential hill. The corresponding Zeeman effect 
on the C >2 internal level states is plotted in the lower part. 
Laser excitations and spontaneous decays are respectively il¬ 
lustrated with red solid line and red dashed line. 


The lifetime of the B-state is 75 ns [371131] and wave¬ 
lengths from its first vibrational level to the X state vi¬ 
brational levels are 541, 598, 667, 753, 863, 1007, 1206 
nm with a corresponding transition strength probability 
of 72, 23, 5, 0.8, 0.1, 2 x 10 -4 , 3 x 10 -5 , 4 x 10 -6 (per¬ 
centage of the FC), from [39, 40 . The lifetime of the 
A-state is 50 ps with wavelengths of 2.53 pm (branching 
ratio of 96%) and 4.57 pm (branching ratio of 4%). In 
order to close the rotational transition cycle, two lasers 
are required for each of these vibrational transitions to 
couple X(N" = 0, 2) to A (TV' = 1, J' = 1/2), see Fig. 
0b). 

For our simulations, the temperature along one axis is 
calculated from the deviation of 50% of the central ve¬ 
locities’ histogram. The so-called 3D-temperature (T^d) 
is the quadratic mean of the three one-dimensional tem¬ 
peratures (x,y,z). 

For our first simulation, we study a possible experi¬ 
ment based on deceleration of an anionic beam. In con¬ 
trast to neutral cold molecules which are difficult to pro- 
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duce at low velocity, and in spite of the development of 
techniques such as velocity filtering, buffer gas cooling 
or decelerators (see list in ESI). an anionic beam can be 
brought to a standstill very easily by an electric field m ■ 
Indeed, a typical beam of C 2 has a current of 1 nA and 
is emitted at 1000 eV with an energy dispersion of 1 eV 
[42]. A 1000 volt potential box can thus decelerate such 
a beam. Furthermore, due to the energy dispersion of 
the beam, the stability of this voltage power supply is 
not critical. Typically 1/100000 of the anions (i.e. 0.01 
pA current) will be decelerated inside the box to below 
an energy of 0.01 meV. This corresponds roughly to 0.1 
K, which is within the capture range of molasses cooling. 

Thus, we propose a very simple deceleration scheme us¬ 
ing a grounded vacuum chamber through which the 1000 
eV beam propagates; this is followed by a chamber at 900 
V with a 3 mm radius hole which focuses the beam (1 
cm downstream) to a smaller hole of R ^ 0.3 mm radius 
of an innermost chamber at 1000 V. The chambers could 
be made of transparent conducting film (such as Indium 
Tin Oxide) or as grids to allow laser cooling. In this fi¬ 
nal chamber the electric field decreases as 0-2 (z/R)~ 3 
along the propagation axis z [43] such that at 1 cm the 
Lorentz force becomes negligible compared to the laser 
cooling force that is typically Kkr /10 ~ 10 -20 N, corre¬ 
sponding to an acceleration of 10 5 m • s -2 . 

We simulate such a 3D molasses cooling with sev¬ 
eral CW lasers, all with 0.5 W power and 2 cm waists 
and wavelengths 10 MHz red detuned from the tran¬ 
sitions (Fig. fTb . We consider here a cooling on the 
X(v" = 0,AT = 2) ++B(v' = 0,7V' * 0) transi¬ 
tion. For the 3D cooling we thus have 6 lasers on 
the X(v" = 0, TV" = 2) ^ B(v' = 0,7V' = 0) tran¬ 
sition, 2 repumping counter-propagating lasers on the 
X(v" = 0, N" = 0) ^B(v' m 0,7V' « 0), plus extra 
repumping lasers: 4 lasers for each higher vibrational 
level of the ground state X up to v" = 4 (these lasers 
are counter-propagating along +/— Z). We start with 
150 particles at 80 mK distributed following a spheri¬ 
cal gaussian distribution (a = 4 mm). The result is 
indicated in Fig. [2ja). The first increase at 1 ms is 
due to the conversion of the Coulomb potential energy 
into kinetic energy. Then, a fast cooling down to 1 mK 
is observed, although many molecules escape from the 
laser waist interaction area. We find that the final tem¬ 
perature and losses depend on the initial density (ini¬ 
tial gaussian radius) because of space charge effects that 
counter the cooling. Note that we do not attempt trap¬ 
ping but only molasses-cool here. Indeed, realizing a trap 
(MOT) would require a magnetic-field gradient produc¬ 
ing a Lorentz force stronger than the laser trapping one, 
especially because the Zeeman effect of B(v' = 0, N' = 0) 
is weak (quadratic behavior). To be feasible, this cooling 
method requires pre-cooling down to few hundreds of mK 
in order to avoid having too fast molecules. 

In the second simulation we study particles in a Paul 
trap. For simplicity, we only consider motion in an har¬ 
monic pseudopotential well K(r) = 1/2 m c j r (x 2 + y 2 ) + 


Doppler Laser cooling simulations 



FIG. 2: (Colour online). Doppler laser cooling results, cycled 
on the Xe^B transition, (a): Deceleration and 3D molasses 
cooling of C 2 ■ Result of the 3D laser cooling simulation of 
an initial bunch of 150 anions at 80 mK: The red dashed 
curve gives the population losses, corresponding to molecules 
which escape from the laser interaction area. Black curve: 
kinetic energy evolution (in temperature unit) of the remain¬ 
ing particles, (b): In black, evolution of the temperature and 
number of C^~ initially at 5K in a Paul trap. The red dashed 
curve gives the population losses, through photo-ionisation or 
decays into high vibrational ground states (v > 4). 


1/2 m cj 2 [44]. We chose 5 K as a typical starting tem¬ 
perature of the beam. This would correspond to having 
carried out a first cooling step, using e.g. Helium as a 
buffer gas. We simulate Doppler cooling on the X f) B 
transitions, with a 600 MHz red detuning, and a spectral 
laser bandwidth of 35 MHz. As previously, we consider 
repumping lasers from the different vibrational levels of 
the ground state X up to X{v" = 4), as shown in Fig. 
[lja). But in this case repumping lasers are only along 
+Z, since particles are trapped. Results are given in Fig. 
|2|b) where cooling down to 60 mK is achieved within 50 
ms. As the photo-detachment cross-section of C</ for the 
B-state is unknown, we use the photo-detachement cross 
section of C 2 , cr ~ 10 -17 cm 2 m as typical value. The 
photo-detachement rate for the B-state is la/hv = 4.3 
s -1 for I = 0.16 W/cm 2 . Within 50 ms of cooling we 
loose only 3% of the molecules by photo-detachement. 
To these photo-detachement losses, we have to add 25% 
of decays to higher vibrational levels of the X-state. The 
losses’ evolution over time is shown in Fig. Igb). Here, 
we would like to emphasize that for both molasses cooling 
and Paul trap simulations, we use a Doppler laser cool¬ 
ing process, with the same lasers. The Paul trap can thus 
serve as a first cooling and can be turned off for further 
cooling using the molasses cooling for low density clouds. 
Photo-detachement and decays in higher vibrational X- 
state levels are thus similar. For the simulations of mo¬ 
lasses cooling however, losses due to motion of molecules 
out of the laser’s area are much more important than 
photo-detachment or decay losses. 

A solution to avoid the losses of through photo¬ 
detachment or decays into high vibrational levels of the 
ground-state is to cool through the A-state. This cooling 
has very similar characteristics (linewidth, wavelength) 
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FIG. 3: (Colour online). Sisyphus cooling curves correspond¬ 
ing to the evolution of the axial temperature and number of 
C 2 initially at 60 K and 100 K in an inhomogenous Penning 
trap and cooled on the XffA transition, with (a) or without 
(b) repumping the X{v" = l)losses. 

to those of La - |2T but with a lighter particle and no 
photo-detachment in the case of C^". 

Our third simulation thus concentrates on cooling and 
trapping in a ^1.5 cm long Penning-like trap using 
Sisyphus-type cooling [28] . The principle is illustrated in 
Fig. Qb): due to the axial motion induced by the electric 
trapping potential (300 qs oscillation time in the simu¬ 
lation) in a Penning-like trap, particles move between 
the high (2 T) and low (0.2 T) magnetic fields at both 
ends of their axial range. A given particle is initially in 
the X(v" = 0,7V 1 2 3 4 5 6 7 8 ' = 0 ,M" = 1/2) state; then, after an 
absorption followed by spontaneous emission, it decays 
towards the X(v" = 0, N" = 0, M" = —1/2) state. More 
than 1 K is removed for each closed absorption-emission 
cycle, as the molecule continuously climbs the magnetic 
potential hill. 

In the last 3D simulation we focus on the axial tem¬ 
perature, set initially at 100 K for a cloud of 200 an¬ 
ions. For our low density plasma there is no coupling 
between radial and axial motions, but the radial shape 
of the plasma reflects the inhomogeneous magnetic field. 
Two cooling lasers at2.5351 qm and 2.5358 qm with 100 
MHz spectral bandwidth each are detuned to be resonant 
at 0.2 T and 2 T, respectively, along the X(v n = 0, TV" = 
0, M" = —l/2)o. 2T ^ A(v f = 0, TV' = 1, J' = l/2)o. 2T 
and X(v" = 0, TV" = 0, M" = 1/2) 2T A(v' = 0, TV' = 
1, J' = 1/2) 2 t transitions. The considered laser power is 
0.05 W for a 1 mm waist. To avoid relying on too many 
lasers, we repump the losses in the X(v" = 0, TV" = 2) 


states with only one single broadband (6000 MHz, 0.05 
W) laser which addresses all the Zeeman-split sub-levels, 
resonant at 0.2 T. 

Results are given in Figj3^a). In tens of ms, the axial 
temperature is cooled down from 60 K to few K. The lost 
population mainly goes to X(y" = 1). As for all simula¬ 
tions, we load the trap using an initial (non thermalized) 
Gaussian velocity distribution. The evolution of this non 
equilibrium system leads to the high frequency veloc¬ 
ity (and thus instantanous temperature) fluctuations in 
both Fig{2|b) and Figj3j We also present an alterna¬ 
tive simulation, of which results are given in Fig. !b). 
Here, decays in the vibrational X(v" = 1) states are re¬ 
pumped: two lasers, at 4.574 pan and 4.593 qm, address 
the X(v" = 1,7V" = 0) A(v' = 0,7V' = 1, J' = 1/2) 

and X(v" = 1,7V" = 2) A(V = 0,7V' = 1, J' = 1/2). 
Both are resonant at 0.2 T, with both spectral bandwidth 
of 6 GHz and power of 0.1 W. 

This simulation requires 2 more lasers but has the ad¬ 
vantage of cooling a greater part of the molecules (less 
than 0.5% of anions fall in the X(y" = 2) levels, within 
80 ms, in comparison to the 60% of losses within 20 ms, 
for the first case without repumpers on X(u" = 1)). 

In conclusion, we have presented several possible decel¬ 
eration and laser cooling schemes for anionic molecules, 
either in free space or trapped, by using Doppler or Sisy¬ 
phus cooling, circumventing the problem of photodetach¬ 
ment. Working with traps can open many possibilities 
due to the long trapping times : it could enable restrict¬ 
ing studies to only ro-vibrational transitions, or working 
with electronic transitions that have long spontaneous 
emission times. Furthermore, laser cooling of molecu¬ 
lar anions followed by laser photo-detachment could be 
used as a source for cold neutral molecules, or as an ideal 
source for electron bunches, since no Coulomb force due 
to ions will be present and ideal uniform elliptical density 
shapes could be realized [45] 06] . 
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